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^ , Abstract 

In this paper, we explore how to constructively manipulate qubits by rotating Bloch spheres. It 
is revealed that three-rotation and one-rotation Hamiltonian controls can be constructed to steer qubits 
when two tunable Hamiltonian controls are available. It is demonstrated in this research that local-wave 
function controls such as Bang-Bang, triangle-function and quadratic function controls can be utilized 
to manipulate quantum states on the Bloch sphere. A new kind of time-energy performance index is 
proposed to trade-off time and energy resource cost, in which control magnitudes are optimized in terms 
of this kind of performance. It is further exemplified that this idea can be generalized to manipulate 
encoded qubits on the controllable subspace. 
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I. Introduction 

Dating from the birth of quantum theory, control of quantum systems is an important issue 
[1]. Quantum control theory has been developed ever since 1980s [2], [3], [4]. Recently, quantum 
information and quantum computation is the focus of reseach[5]. A great progress has been made 
in the domain of quantum control[6], [7], in which the controllability of quantum systems is a 
fundamental issue. The different notations of controllability have been explored in [8], [9], [10], 
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[11], [12], [13], [14]. Specially, the controllability of quantum open systems has been studied 
by some researchers [15], [16], [17], [18]. It is quite well known that quantum open systems are 
not open-loop controllable but there may exist decoherence-free subsystems or subspaces[19], 
[20], [21], [22], [23]. The works on encoded universality [24], [25] further enhance the belief that 
one can manipulate quantum information on the encoded subspace. Optimal control theory has 
also been successfully applied to the design of open-loop coherent control strategies in physical 
chemistry [26], [27], [28]. Recently, time-optimal control problems for spin systems have been 
solved to achieve specified control objectives in minimum time[29], [30], [31]. On the other hand, 
the challenge of open-loop control is to design external fields or potentials acting as model-relied 
controls. The main strategies for open-loop control design seem to be based either on geometric 
ideas or more formally Lie group decompositions, as in[32], [33], [34], [35], [36]. 

In this paper, we explore how to constructively manipulate qubits or encoded qubits based on 
the geometric parametrization of qubits when two tunable Hamiltonian controls are available. It is 
demonstrated that one can not only design 3-rotation Hamiltonian controls to manipulate qubits, 
but can also construct 1 -rotation Hamiltonian controls to steer qubits by carefully choosing a 
rotation axis. It should be underlined that local wave controls can be constructed to manipulate 
qubits corresponding to each rotation. Furthermore, we proposed a new kind of time-energy 
performance index 



where E(u(t)) is the energy cost of control at time t, tf is free terminal time, and A is introduced 
as a ratio parameter to trade-off the cost of time and energy resource. It has also been discussed 
in [36] how to optimize 3-rotation Bang-Bang controls to transfer quantum state in terms of this 
kind of time-energy performance. In this paper, we comprehensively discuss how to optimize 
control magnitudes in terms of this kind of time-energy performance for both 3-rotation and 
1-rotation controls, and present optimal Bang-Bang, triangle-function and quadratic function 
controls, respectively. 

The rest of this paper are organized as follows. In Sect. II, we present prerequisite for further 
discussion. It is illustrated in Sect. Ill how to manipulate qubit by 3-rotation Bang-Bang, triangle- 
function, and quadratic function controls. The optimal controls are further presented in the sense 
of time-energy performance. It is also revealed in Sect. IV that one can utilize three kinds of 
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local-wave controls to manipulate qubits just by one-times rotation. The paper concludes with 
Sect. V. 

II. Prerequisite 
Consider a controlled qubit governed by the equation 

j t \m) = ~H(t)\m) = ~Mt)cr z + Uy (t)a y ]\m) w 

where a z = I 2 — 2|1)(1| and a y = i\l)(0\ — i\0)(l\. For simplicity, we set H = 1. 

Denote \u + ) = cos^lO) -Msin^|l); |it_) = sin^|0) — icos^|l). It is interesting to point 
out that if H(t) = f(t)[cos6 u a z + sin 9 u a y ], then one can express the Hamiltonian H(t) as 
H(t) = f(t)a z where cr" = \u + )(u + \ — 




As shown in FigfH one can not only choose 3-rotation control functions to steer the controlled 
qubit system from an arbitrary initial state to another arbitrary target state, but can also construct 
1-rotation control function f(t) to achieve the same goal. 

In this paper, we will just concentrate on three kind of local wave-functions: a piece-wise 
constant function (Bang-Bang control), a triangle-function and a quadratic function. 

Denote the triangle function Ut(£; to, h, L) and the quadratic function «q(£; t , ti, L) respec- 
tively as follows: 

ur(*;tb,ti,£) = | -j^-it-h) fe[*±Mi) (3) 

otherwise 
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and 

u Q (t;t ,t u L) = \ ^ 2 [01J (4) 

L otherwise 

where both -u r (t;t ,£i, L) and UQ(t;t ,ti : L) are nonzero only when t G (to,£i)> and take the 
maximum magnitude L at time t -^ LL . It should be underlined that the pulse area of the control 
pulses is the key control variable for geometric control and the pulse area inequality for Bang- 
Bang, triangle-function and quadratic function controls is given as 

rti rt\ rti 

/ UT{t]tQ,ti,L)dt < / UQ(t;t ,ti, L)dt < / Ldt = L(ti — t ) (5) 

J to J t() J to 

Furthermore, it is worth pointing out that f* 1 L 2 dt = L 2 {t\ — t ) and 

£(u T (*;*o, *!,£)) = / |ur(*;*o,<i,^)| 2 ^ = ^ 2 (*i-*o) (6) 

Jto 6 

and 

E{u Q {t-t G ,t u L))= I \u Q (t;t ,t u L)\ 2 dt = A^-to) ( 7 ) 

Jto 15 

Remark: We would like to further emphasize that that one can construct both 3-rotation and 
1-rotation local wave-function controls to manipulate qubits if H{t) = u x (t)a x + u y (t)a y or 
H{t) = u z (t)a z + u x (t)a x . In other words, 3-rotation and 1-rotation controls can be constructed 
as long as two tunable Hamiltonian controls are available. 

III. Manipulate qubit by three rotations 
Consider a controlled qubit governed by the equation 

j t m)) = -i(u z (t)a z + u y (t)a y )\i)(t)) (8) 

with an initial state \i{) ) = cos^lO) + e^°sin^|l) and a target state \ip s ) = cos^|0) + 
e i<f> s s j n For the sake of the following analysis in this section, we denote (j) 0m = min(0 o , 2n— 
4>o), 4> sm = min(0 s , 2n - 4> s ), 9 0s = \6 S - 6 \ and = (f> 0m + %s + <j>sm- 

In this section, our control goal is to find tj and some form of controls {u z (t),u y (t) : < 
t < tj} so that 

m f )) = \ip s ) (9) 

by three rotations about z— axis, y— axis and z— axis, respectively. Furthermore, we hope to 
optimize control magnitude in terms of the performance CQ) where E(u(t)) = \u z (t)\ 2 + \u y (t)\ 2 
and A > 0. 
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A. 3-rotation Bang-Bang controls 

In this subsection, we will discuss how to manipulate quantum system © by three-rotation 
Bang-Bang control. According to the properties of Pauli matrices[5], we choose the piecewise 
constant controls {u z (t),u y (t) : < t < tf\ as follows: 

sign{(f) - tt)M z1 te[0,t{) 

u*(t) = j o te[t u t 2 ) (10) 

sign(7r - <p s )M z2 t G [t 2 , t f ] 



and 



o te[o,ti) 
u y (t) = | S ign(9 s - 9 )M y t G [t 1} h) (11) 
te[h,t f ] 



where h = t 2 = ^- + t, and t f = ^+t 2 . 

After some calculations, we have \ip(ti)) = cos ^|0)+sin \ip(t 2 )) = cos ^|0)+sin 
and \ip(t f )) = cosf |0> + sin f |1). 

Next, our task is to choose M z \, M z2 and M y to minimize the performance (0Q). It can be 
demonstrated that 

J = Mfe + jgft + fe) + + ^ + > VX^e d2) 

where the equality holds only if M z \ = M z2 = M y = y/\. 

If only bounded Bang-Bang controls with bound Lb are permitted, then the optimal controls 
are given as: 

sign(<p - n)L* B £g[0,£J) 

<(*) = { o te[n,t*) (13) 

sz#n(7T - (f) s )L* B t G \t* 2 ,t* f ] 

and 

o *e[0,*i) 

= { sign(9 s - Q )L* B t G [t\,t* 2 ) (14) 

* G [t 2 ,t* f ] 

where tj = fg|, i* = ^ + t*, £} B = |g| + 1* = and L£ = min( v / A, L B ). Furthermore, 

the corresponding optimal performance is J B = (^rr + It is interesting to underline 
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that J* B = A ■ t* fB + E* B where E* B = |L^E^, and t * fB ■ E* B = only depends on the location 
of both initial and target states on the Bloch sphere. 

If unbounded Bang-Bang controls are permitted, we have L* B = V% t* fB = ^4=, E* B = 
and J* B = therefore t* fB ■ E* B = ±£^. 

B. 3-rotation triangle-function controls 

In this subsection, we will first demonstrate that the target state \ip(tf)) = \ip s ) can be achieved 
from the initial state \ij) ) by the following three-rotation triangle-function controls: 

u z (t) = sign((j)Q - ir)u T (t; 0, h, M zX ) + sign(n - 4> s )u T {t] t 2 , t f , M z2 ) (15) 

and 

u y (t) = sign(9 s - 9 )u T (t; h, t 2 , M y ) (16) 

where t x = |gt, t 2 = f£ + h and t f = + t 2 . It can be proved that = cos &|0) + 

sinf |1), \iP(t 2 )} = cosf |0)+sinf |1), and |^(t/)) = cos f |0) + sin 

Subsequently, our task is to select magnitude M z ±, M z2 and M y to minimize the performance 
(OQ). It can be demonstrated that 

^~Mzi 'My ~M^2' \ 3 + ^~~ + 3 > — "TS" ^ K ' 

where the equality holds only if M z \ = M z2 = M y = \/3X. If only bounded triangle-function 
controls with bound L B are permitted, then the optimal 3-rotation triangle-function controls are 
given as: 

u*(t) = sign(4>o — ir)uT{t; 0, t*, L^) + sign^ir — <fi s )uT{t; t* 2 , t% T , Ly) (18) 

and 

u* y {t) = sign{6 s - 9 )u T (t; **, L* T ) (19) 

where t\ = t* = ^ + t*, t* fT = 5^ and L* T = mm(y/3X, L B ). Furthermore, the optimal 
performance corresponding to bounded triangle-function control is J£ — (^r + -y : )E^. It is 
interesting to underline that J£ = A ■ t * fT + Ej, with E? = ^L^E^g, and t* fT ■ E^ = only 
depends on the location of both initial and target states on the Bloch sphere. 

If unbounded triangle-function controls are permitted, then we have t*j T = -^=, U T = V3A, 
Ej, = : ^Y l( pg and J|> = ^-^E^, thus tj T ■ E^ = |E^. 
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C. 3-rotation quadratic function controls 

In this subsection, it is demonstrated that the target state \ip(tf)} = \ip s ) can be achieved from 
the initial state \ip Q ) by the following quadratic controls: 

u z {t) = sign((f>o - n)u Q (t\ 0, t l7 M zX ) + sign(ir - 0s)u Q (t; t 2 ,t f , M z2 ) (20) 

and 

u y (t) = sign(6 s - 9 )u Q (t; t u t 2 , M y ) (21) 

where h = fjjg^, t 2 = ff^ + h and t f = ffg + t 2 . It can be confirmed that = 
cosf |0) + sinf |1), \ifj(t 2 )) = cosf |0)+sinf |1), and |^(*/)) = cos ||0) + e 1 ^ sin 

Next, our task is to choose magnitude M z \, M z2 and M y to minimize the performance ©• 

After some calculations, we have 

7 \ / 3<ftp m , 36q 3 , 34> S m \ I /2Mzi0om , 2My9o 3 . 2M z2 4>sm \ \ V30A y ('99 N > 

J - A ^4M^ + 4M^ + 4M^J + I 5 + + 5 J - "H -2 ^ ^ 

where the equality holds only if M zl = M z2 = M y = ^f^. 

If only bounded quadratic function controls with bound L B are permitted, the optimal 3- 
rotation bounded quadratic-function controls are given as: 

u* z (t) = sign((j) - rr)u Q (t; 0, t* v L* Q ) + sign(n - <f> s )u Q (t; t*, t* f , L* Q ) (23) 

and 

u* y (t) = sign(9 s - o )u Q (t;tl,t* 2 ,L* Q ) (24) 
where t\ = t* 2 = fg£ + t ls t* fQ = g# and L* Q = min(^p, L B ). Moreover, the optimal 

Q Q Q 2Z/* 

performance corresponding to bounded control is Jq = (^4- -\ — g 2 -)!^. It is interesting to 
underline that J Q = A ■ t* fQ + E Q where E Q = |Lq£^, and tj Q ■ Eg = only depends on 

the location of both initial and target states on the Bloch sphere. 

If unbounded quadratic function controls are permitted, we have t*j Q = ^=|, L*q = -^-f^, 
E* Q = ±§*£^ and J* Q = ^p£^, therefore t* fQ ■ E* Q = 

Remark: 1. When unbounded controls are permitted, it has been demonstrated in this section 
that J* B < J Q < J?, E* B < E* Q < E?p and t* fB < t* fQ < t* fT , therefore we have t* fB -E* B < 

+* TP* 7-7 1 * 
fQ Q < l fT'-^T- 

2. Even when only bounded controls are permitted, the above inequalities are valid for all A 
and L B except that E* B < E Q < E^ does not hold for some A and L B . 
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IV. Manipulate qubits just by one rotation 

Reconsider the controlled qubit ([8]) with both the same initial and target states given in the 
Section III. In this section, our control goal is to find tf and some form of controls {u z (t), u y (t) : 
< t < tf} so that \ip(tf)} = \ip s ) is attained just by one rotation. Furthermore, we hope to 
choose {u z (t),u y (t) : < t < tf} to minimize the performance CO)- 

Choose H(t) = f(t) (cos 9 u a z + sin 9 u a y ) with 9 U E [0,7r] so that the following equation holds 

sin 9 U sin 8q sin (fio + cos 9 U cos 9q = sin 9 U sin 9 S sin <p s + cos 9 U cos 9 S (25) 

Since |0) = cos^-|m + ) + ism.^f\u-)\ |1) = — ism.^f\u + ) + i cos ^|it_), the initial and target 
states can be expressed in terms of the new basis \u + ) and \uJ) as follows 

|^) = cosS|u+) + e < ^sin^|u_) (26) 



and 



where 



and 



and 



= cos^K) + e^sin^|«_) (27) 



cos -f- = \ - + - [sin 9 U sin 9 sin (f) + cos 9 U cos 9q] (28) 



2 V 2 2' 



$ = -Z(cos f cos \ - ie irt>0 sin f sin ^) 
+Z(i cos f sin \ + ie**° sin f cos \) ± 2n 7r 



0f = -Z(cos ^ cos ^ - ie i<t,a sin ^ sin If) 
+Z(i cos ^ sin If + ze i0a sin ^ cos ^) ± 2n s 7r 



(29) 



(30) 



It is easy to prove that one can choose the suitable integers n and n s so that 4>q ,4>f E [0, 2tt). 

Remark: 1. We would like to point out that the initial and target states have the same angle 
6f Q about the control Hamiltonian axis as shown in FigQ] 

2. For the sake of the analysis, we introduce (pf = min(|0^ — 0^|,27r — \<j)f — 0^|) with 
4>^ E [0,7r). It should be underlined that (p^ depends not only on the location of both initial 
and target states on the Bloch sphere, but also on the Hamiltonian H(t), i.e., the y — z plane. 
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A. 1-rotation Bang-Bang controls 

In this subsection, we will discuss how to manipulate the quantum system ([8]) by Bang-Bang 
control. According to the aforementioned analysis in the section, we can choose the piecewise 
constant controls {/(£) : < t < tf} as follows: 

( M ub ,te[0,t f ) i/O±2fc7r<(0f-^)<7T±2fc7r (3i) 
I -M ub , t G [0, t f ) ifn ± 2kn < (0f - 0^) < 2n ± 2/cvr 

where */ = mt b - 

Subsequently, our task is to choose M u b to minimize the performance ([I]) where E(u(t)) = 
\f(t)\ 2 and A > 0. 
After some careful calculations, we have 



J = + ^ > v^ (32) 

where the equality holds only if M u t, = \f\. 

If only bounded Bang-Bang controls with bound Lb are permitted, then the optimal controls 
are given as: 

/ cos9 u L B ,te[0,t* fB ) if0±2kn<(^-^)<n±2kn 
u z (t) = < (33) 

I - cos 9 U L* B , t G [0, t* fB ) if Ti ± 2krr < (cf)f - 0^) < 2tt ± 2A;tt 

and 

/ sin^,te[0 ; ^ B ) i/O±2/ OT <(0f-0 ( f)<7T±2/ OT 
u At) = \ (34) 

I - sin 6 U L* B , t G [0, t* fB ) if n ± 2A;tt < (0f - 0f) < 2tt ± 2/br 
where tf R = ttttt and L* R = mm(\/A, ; — 4r 5—). The optimal performance corresponding 

J 11j b max(cos y,siii -^-) ' *- 1 

to bounded Bang-Bang controls is J B = {^JF + ^-) ( t ) %- It ls interesting to emphasize that optimal 
performance can be expressed as J* B = A ■ t* fB + E* B with E* B = 4f</^, and E* B ■ t* fB = \((f)^ ) 2 
where (p^ is independent of A. 

If unbounded Bang-Bang controls are permitted, then L* B = V% t* fB = E* B = ^<pf , 
and J* B = V\<J>h. Therefore, E* B ■ t) B = ±(0f o ) 2 . 
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B. 1-rotation triangle-function controls 

In this subsection, we will explore how to construct one-rotation triangle-function controls 
{f{t) '■ < t < tf} to achieve the target state from the initial state. We can select 



u T (t; 0,t f , M ut ) ifO ± 2kii < (0f - <$) < vr ± 2k 



7T 



f(t) = t w ' J1 ' " ~ KTS TU/ (35) 

I -u T (t; 0, t f , M ut ) if TV ± 2kn < (0f - <f>%) < 2n ± 2A;tt 

where = jf-. In other words, {w 2 (£), u y (t)} can be constructed as follows: 

f ^(t;0,|& M^cosfl* i/O±2fc7r<(0f-0f)<7r±2A;7r 
w) = i Th _ _ . ..„ (36) 

and 



-mt(*; 0, M„t) cos 9 U if -K ± 2A;tt < (0f - 0*) < 2tt ± 2A;tt 

f u T (t;0,^,M ut )cosfl u i/0 ± 2A;tt < (0f - 0^) < tt ± 2&tt 
I -u T {t;0,^L,M ut )cos9 u ifn ± 2A;tt < (0f - 0f) < 2tt ± 2A;vr 
Next, our task is to optimize magnitude M ut in terms of the performance ([I]). It is easy to 
demonstrate that 

J = a€ + > ^ < 38 > 



where the equality holds only if M Mt = y3A. 

If only bounded triangle-function controls with bound L B are permitted, then the optimal 
controls are given as: 

f u T {t;0,t* fT ,L* T )cos6 u ifO ± 2fc7r < (0f - 0f) < tt ± 2kir 

u z (t) = < (39) 

I -u T (t;0,t* fT ,Ur)cos9 u ifn±2kn<((j)^-(f)^)<27i±2kTT 

and 

= / u T(t;0,t} T ,L* T )sm6 u if 0±2kw <n±2ter ^ 

V I -u T (t; 0, £} T , Ly) sin 6 U if tt ± 2/ctt < (0f - 0f ) < 2tt ± 2A;tt 

where £y T = and L* T = min(v / 3A, ma ( cos L ^ S i n ^ j ) • The °pti ma l performance corresponding 
to bounded control is = (p- + -^)4>%. It is interesting to emphasize that optimal performance 
corresponding to bounded triangle-function controls can be expressed as J£ = A ■ t* fT + Ej, with 

E?p = ^L^(f)f , and E? ■ t* fT = §(0^) 2 where (f)f is independent of A. 

<t> H i — 

If unbounded triangle-function controls are permitted, then t*j T = -^=, L\ = v3A, E^ = 
^p0f o and J* T = 0f o . Therefore E* T • t* fT = §(0f o ) 2 . 
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C. 1-rotation quadratic-function controls 

In this subsection, we will explore how to construct quadratic controls {/(£) : < t < tf} to 
achieve the target state from the initial state. We can choose 

f(t) = { U ^ W /' M ->) if0±2for<(,<ff-<ffi)<ir±2for ^ 
I -u Q (t;0,t f ,M uq ) ifn ± 2kn < (</>f - <ftf) < 2tt ± 2kn 

where tf = j^f-- In other words, the quadratic controls {u z (t),u y (t)} are given as follows: 

U (t) = J ^°>l^ M ^) COS ^ ± 2fc7T < (0f - O H ) < 7T ± 2fc7T ^ 

I -ug(t; 0, M nq ) cos ft u ± 2fc7r < (0f - <) < 2tt ± 2A;tt 

and 

( MQ (t;0,^,M n9 )sm^ z/0 ± 2kn < (0f - <ftf) < tt ± 2&7T 

I -u (f; 0, |g, M U9 ) sin0 u */tt ± 2A;tt < (0? - 0*) < 2tt ± 2A;tt 

Next, our task is to choose magnitude M ug to minimize the performance (OQ). After some 
calculations, we further obtain 

j = X ^L + > VfE^H (44) 

where the equality holds only if M uq = 

If only bounded quadratic controls with bound Lb are permitted, then the optimal controls 
are given as: 

„*(£) = / «q(*;M/q,£q) cos U i/O±2fc7r<(0f-^ r )<7r±2fc7r (45) 

I -u Q (t;0,t* fQ ,L* Q )cos9 u ifn ± 2kn < (0f - < 2tt ± 2/ett 

and 

«*(*) = / UQ(t;O,t} Q ,L£)sin0 u 2/0 ± 2/ctt < (0f - O H ) < tt ± 2A;tt ^ 

" \ -u Q (t;0,t fQ , L* Q ) sin 9 U ifn ± 2/br < (0f - 0f) < 2tt ± 2A;tt 

where t} Q = and Lq = min(^p, max(c J^ sin ^ } )- The optimal performance corresponding 
to bounded control is Jq = (^- H — g^X^)- It is interesting to emphasize that optimal 
performance corresponding to unbounded quadratic controls can be expressed as J Q = X-^q+Eq 
with E* Q = lL* Q( pu, and E* Q ■ t* fQ = ^(<^) 2 where (j)f is independent of A. 

If unbounded quadratic controls are permitted, then t* fQ = -j=, L* Q = ^f^, E Q = ^p</>^, 
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D. Further discussions 

1. When unbounded controls are permitted, we have J* B < J Q < Jy, E* B < E Q < E? and 
t*f B < t* fQ < t*f T , therefore we have t* fB -E B < t* fQ -E Q < t* fT -E^. 

2. Even when only bounded controls are permitted, the aforementioned inequalities are valid 
for all A and L B except that the inequality E* B < Eq < E? is invalid for some A and L B . 

3. When one fixed Hamiltonian and another tunable control Hamiltonian are available, only 
1— rotation Bang-Bang control can be designed to transfer the qubit from the initial state to 
the target state. For example, if Hit) = (a z + u y (t)a y )\tp(t)) and sin^osin^o ^ sin 9 S sin (f> s , 
one may be able to construct 1— rotation Bang-Bang control to achieve the target state. When 
unbounded Bang-Bang control are available, one should choose u y (t) = tan# u where tan# u = 

cose a -cose _ When only bounded Bang-Bang controls with the bound L R are available, 
1— rotation bounded Bang-Bang control can be constructed only if L B > \ tan^ u |. This result is 
in interesting contrast with the recent research[39]. 



V. Discussions and conclusions 

At first, we would like to point out that the three-rotation and one-rotation control design 
methods can be generalized to manipulate encoded qubit on controllable subspace of both closed 
and open quantum systems. 

For example, let us consider a controlled 2-qubit system which is governed by the equation 

^\m) = -{H{u{t))\m) (47) 

where H(u(t)) = u zih (t)ai 1] ' +u hZ2 (t)Ip ] ®af ] +u yiX2 {t)a y l) ®af ] +u Xiy ^t)a { i ] ®a y 2) . 
Under the above condition, an encoded qubit basis can be given as { 1 L ) = |0 1 l 2 ),|li / ) = |li0 2 )}. 
Denote the encoded subspace, which can be expanded by the encoded state basis {\0l)i as 
E L . It is interesting to underline that for any pure state \iPe)^E l , one can obtain its geometric 
parametrization in terms of {|0l) = IO1I2) and = | I1O2)}. Denote a\ = |0l)(0l| — 
\1 L )(1 L \ = KaP^I^-lP^aP) and a L y = i|li)<0 L |-i|0 £ )(l £ | = \{a y l) ®a { x 2) -a { } ] ®af ] ) 
By setting u Zl j 2 (t) = —u IlZ2 (t) = \u^(t) and u yiX2 (t) = —u xiy2 (t) = ~Uy{t), one can express 

j t \m) = ~(u L z {t)o- L z + u L y {t)o L y )\m) (48) 



the equation PTI) as 
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For an open quantum system, its dynamics equation is in general rather difficult to gain. 
However, in many practical situation, quantum dynamical semi-group master equation[37], [38] 
is an appropriate way to describe the evolution of the quantum open system as follows 

^ = ~[H(<t)),P] + Hp) (49) 

where Lindbladian is: 

1 - 

L (P) = 2 E a ^ P F ft + &P> F N (50) 

and H(u(t)) is the system Hamiltonian, the operators Fi constitute a basis for the TV-dimensional 
space of all bounded operators acting on H, and are the elements of a positive semi-definite 
Hermitian matrix. 

If H(u(t)) = u Zlh (t)a[ 1] <g> 4 2) + u hZ2 (t)lP ® of + u yiX2 {t)a y l) ® a x 2) + u Xiy2 {t)o { x ] <g> a y 2) 
and L(|'0b)(V'b|) = for any pure state \iPe)^E l , then, for p = \iPe){^e\ with \^e)^E l , 
Eq.(|49l) is further reduced to Eq. (g8]) because L(|-0k) = 0. 

So far, it has been demonstrated in this research that one can utilize various local wave-function 
controls including Bang-Bang controls, triangle-function controls and quadratic-function controls 
to manipulate qubits and encoded qubits on controllable subspaces for both open quantum 
dynamical systems and uncontrollable closed quantum dynamical systems when two tunable 
Hamiltonian controls are available. Furthermore, we discuss how to design control magnitude in 
terms of a kind of time-energy performance. It is demonstrated that optimal Bang-Bang controls 
have the best performance and optimal triangle-function controls have the worst performance 
among three kinds of control schemes. It is the pulse area inequality for three controls given in 
Eq. © who makes the performance difference. It should be emphasized that one can introduce a 
ratio parameter A to trade-off between time and energy resource cost, but the product of time and 
energy cost is an invariance under different A for each kind of controls due to the characteristic 
of geometric control. 

It is well known that low-capacitance Josephson tunneling junctions offer a promising way to 
realize qubits for quantum information processing [40] and two tunable Hamiltonian controls are 
available in this application. Therefore this research implies that one can constructively adjust 
gate voltages or magnetic fields to manipulate qubits based on either charge or phase (flux) 
degrees of freedom . 



June 17, 2011 



DRAFT 



IEEE TRANSITION ON AUTOMATICAL CONTROL, VOL. X, NO. XX, XXXX 20XX 



14 



References 

[1] A. Blaquiere, S. Diner and G. Lochak, (edit) Information Complexity and Control in Quantum Physics, Springer- Verlag, 
New York, (1987) 

[2] G. M. Huang, T. J. Tarn and J. W. Clark, "On the Controllability of Quantum Mechanical Systems", J. Math. Phys., vol.24, 
pp. 2608, (1983) 

[3] C. K. Ong, G. Huang, T. J. Tarn and J. W. Clark, "Invertibility of Quantum-Mechanical Control Systems", Math. Sys. 
Theor. vol.17, pp. 335, (1984) 

[4] J. W. Clark, C. Ong, T. J. Tarn and G. M. Huang, "Quantum non-Demolition Filters", Math. Sys. Theor. vol.18, pp. 33, 
(1985) 

[5] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge, Cambridge University 
Press, (2000) 

[6] D'Alessandro, Domenico Introduction to quantum control and dynamics, CRC Press, (2007) 

[7] D. Dong, and I. R. Petersen,"Quantum control theory and applications: a survey", IET Control Theory AppL, 4(12), pp. 
2651-2671 (2010) 

[8] V. Ramakrishna, M. V. Salapaka, M. Dahleh, H. Rabitz and A. Peirce, "Controllability of molecular systems", Phys. Rev. 
A, 51, 960 (1995) 

[9] R. B. Wu, T. J. Tarn and C. W. Li, "Smooth controllability of infinite dimensional quantum mechanical systems", Phys. 
Rev. A, 73, 012719 (2006) 

[10] S. G. Schirmer, H. Fu and A. I. Solomon, "Complete controllability of quantum systems", Phys. Rev. A, 63, 063410 (2001). 
[11] C. B. Zhang, D.Y. Dong and Z.H. Chen, "Control of noncontrollable quantum systems: a quantum control algorithm based 

on Grover iteration", J. Opt. B: Quantum Semiclassical Opt., 7, S313-S317 (2005) 
[12] F. Albertini and D. D'Alessandro, "Notions of controllability for bilinear multilevel quantum systems" IEEE Trans. Autom. 

control, 48(8), 1399 (2003) 

[13] R. Wu, A. Pechen, C. Brif and H. Rabitz "ontrollability of open quantum systems with Kraus-map dynamics", J. Phys. 

A:Math. Theor, 40 5681-5693 (2007) 
[14] G. Turinici and H. Rabitz, "Quantum wavefunction controllability", Chem. Phys. 267(1-3), 1-9 (2001) 
[15] S. G. Schirmer, I. C. H. Pullen and A. I. Solomon, controllability of quantum systems, in Hamiltonian and Lagrangian 

Methods in Nonlinear Control, Proceedings of the second IFAC Workshop, Seville, Spain, 2003, edited by A. Astolfi and 

A. J. van der Schaft (Elsevier Science Ltd., New York, 2003), pp. 311-316 
[16] C. Altafini, "Coherent control of open quantum dynamical systems", Phys. Rev. A, 70, 062321 (2004) 
[17] S. Lloyd and L. Viola,"Engineering quantum dynamics", Phys. Rev. A 65, 010101(R) (2001) 

[18] M. Zhang, H.-Y. Dai, X. C. Zhu, X. W. Li and D. Hu, "Control of the quantum open system via quantum generalized 

measurement", Phys. Rev. A 73, 032101 (2006) 
[19] P. Zanardi and M. Rasetti,"Noiseless Quantum Codes", Phys. Rev. Lett. 79, 3306 (1997) 

[20] D. A. Lidar, I. L. Chuang and K. B. Whaley, "Decoherence-Free Subspaces for Quantum Computation", Phys. Rev. Lett. 
81, 2594 (1998) 

[21] J. Kempe, D. Bacon, D. A. Lidar and K. B. Whaley, "Theory of decoherence-free fault-tolerant universal quantum 

computation", Phys. Rev. A, 63, 042307 (2001) 
[22] M. Zhang, H.-Y. Dai, G.-H. Dong, H.-W. Xie and D. Hu, "Controllable subsystems of quantum dynamical systems", 

Quantum Information and Computation, 1, 469C478 (2007) 



June 17, 2011 



DRAFT 



IEEE TRANSITION ON AUTOMATICAL CONTROL, VOL. X, NO. XX, XXXX 20XX 



15 



[23] F. Ticozzi and L. Viola, "Quantum Markovian subsystems: invariance, attractivity and control", IEEE Transaction on 

Automatic Control, 53(9), 2048-2063 (2008) 
[24] J. Kempe, D. Bacon, D. P. DiVincenzo and K. B. Whaley, "Encoded universality from a single physical interaction", 

Quantum Information and Computation, 1, 33-55 (2001) 
[25] J. Vala and K. B. Whaley, "Encoded universality for generalized anisotropic exchange Hamiltonians", Phys. Rev. A, 66(2), 

022304 (2002) 

[26] M. Shapiro and P. Brumer, Principles of the quantum control of molecular processes, John Wiley and Sons, Inc. (2003) 
[27] A. P. Peirce, M. Dahlen and H. Rabitz, "Optimal control of quantum-mechanical systems: existence, numerical 

approximation, and applications", Phys. Rev. A, 37, 4950-4964 (1988) 
[28] M. Dahlen, A. P. Peirce and H. Rabitz, "Optimal control of uncertain quantum systems", Phys. Rev. A, 42, 1065-1079 

(1990) 

[29] N. Khaneja, R. Brockett and S. J. Glaser, "Time optimal control in spin systems", Phys. Rev. A, 63, 032308 (2001) 
[30] U. Boscain, G. Chariot, J. P. Gauthier, S. Guerin and H. R. Jauslin, "Optimal control in laser-induced population transfer 

for two- and three-level quantum systems", J. Math. Phys., 43, 2107-2132 (2002) 
[31] U. Boscain and P. Mason, "Time minimal trajectories for a spin 1/2 particle in a magnetic field", J. Math. Phys., 47, 

062101 (2006) 

[32] M. Reck, A. Zeilinger, H. J. Bernstein and P. Bertani, "Experimental realization of any discrete unitary operator", Phys. 
Rev. Lett., 73, 58C61 (1994) 

[33] V. Ramakrishna, K. L. Flores, H. Rabitz and R. j. Ober, "Control of a coupled two-spin system without hard pulses" Phys. 
Rev. A 62, 053409 (2000) 

[34] S. G. Schirmer, A. D. Greentree, V. Ramakrishna and H. Rabitz, "Constructive control of quantum systems using 

factorization of unitary operators", J. Phys. A: Math. Gen., 35, 8315C8339 (2002) 
[35] R. Cabrera, T. Strohecker and H. Rabitz, "The canonical coset decomposition of unitary matrices through Householder 

transformations", J. Math. Phys., 51, 082101 (2010) 
[36] W. Zhou, S. G. Schirmer, M. Zhang and H.-Y. Dai, "Bang-bang control design for quantum state transfer based on 

hyperspherical coordinates and optimal time-energy control", J. Phys. A: Math. Theor., 44, 105303 (2011) 
[37] G. Lindblad, "On The Generators Of Quantum Dynamical Semigroups", Commun. Math. Phys., 48, 119 (1976). 
[38] R. Alicki and K. Lendi, Quantum Dynamical Semigroups and Applications, Lecture Notes in Physics Vol. 286, Springer- 

Verlag, Berlin (1987) 

[39] D. Dong, L. James and I. R. Petersen, "Robust incoherent control of qubit systems via switching and optimisation", Int. 
J. Control, 83, 206-217 (2010) 

[40] Y. Makhlin, G. Schon and A. Shnirman, "Quantum-state engineering with Josephson-junction devices", Rev. Mod. Phys. 
73, 357 (2001) 



June 17, 2011 



DRAFT 



